Identification, Expression and Evolutional Analysis of Two cyp19-like Genes in Amphioxus.

The mechanism of sex determination and differentiation in animals remains a central focus of reproductive and developmental biology research, and the regulation of sex differentiation in amphioxus remains poorly understood. Cytochrome P450 Family 19 Subfamily A member 1 (CYP19A1) is a crucial sex differentiation gene that catalyzes the conversion of androgens into estrogens. In this study, we identified two aromatase-like genes in amphioxus: cyp19-like1 and cyp19-like2. The cyp19-like1 is more primitive and may represent the ancestral form of cyp19 in zebrafish and other vertebrates, while the cyp19-like2 is likely the result of gene duplication within amphioxus. To gain further insights into the expression level of these two aromatase-like, we examined their expression in different tissues and during different stages of gonad development. While the expression level of the two genes differs in tissues, both are highly expressed in the gonad primordium and are primarily localized to microsomal membrane systems. However, as development proceeds, their expression level decreases significantly. This study enhances our understanding of sex differentiation mechanisms in amphioxus and provides valuable insights into the formation and evolution of sex determination mechanisms in vertebrates.

A novel C-type lectin protein (BjCTL5) interacts with apoptosis stimulating proteins of p53 (ASPP) to activate NF-κB signaling pathway in primitive chordate.

C-type lectin proteins (CTLs), a group of pattern recognition receptors (PRRs), play pivotal roles in immune responses. However, the signal transduction and regulation of CTLs in cephalochordates have yet to be explored. In this study, we examined the composition of CTLs in Branchiostoma japonicum, identifying a total of 272 CTLs. These CTLs underwent further analysis concerning domain arrangement, tandem and segmental duplication events. A multidomain C-type lectin gene, designated as BjCTL5, encompassing CLECT, KR, CUB, MAM, and SR domains, was the focal point of our investigation. BjCTL5 exhibits ubiquitous expression across all detected tissues and is responsive to stimulation by LPS, mannose, and poly (I:C). The recombinant protein of BjCTL5 can bind to Escherichia coli and Staphylococcus aureus, inducing their agglutination and inhibiting the proliferation of S. aureus. Yeast two-hybrid, CoIP, and confocal immunofluorescence experiments revealed the interaction between BjCTL5 and apoptosis-stimulating proteins of p53, BjASPP. Intriguingly, BjCTL5 was observed to induce the luciferase activity of the NF-κB promoter in HEK293T cells. These results suggested a potential interaction between BjCTL5 and BjASPP, implicating that they involve in the activation of the NF-κB signaling pathway, which provides an evolutionary viewpoint on NF-κB signaling pathway in primitive chordate.

Molecular and functional characterization of a type-1 cystatin in amphioxus (Branchiostoma japonicum).

Cystatins comprise a vast superfamily of evolutionary conserved proteins, predominantly recognized for their roles as endogenous inhibitors by regulating the activity of cysteine proteases. Emerging lines of research evidence also provides insight into their alternative roles in a spectrum of biological and pathological processes, including neurodegenerative disorders, tumor progression, inflammatory diseases, and immune response. Nowadays, various type-1 cystatins (stefins) have been demonstrated among a variety of discovered vertebrate groups, while little is known about the related homologue in cephalochordate amphioxus, which are repositioned at the base of the chordate phylum. In the present study, a single type-1 cystatin homologue in Branchiostoma japonicum was first successfully cloned and designated as Bjcystatin-1. The deduced Bjcystatin-1 protein is structurally characterized by the presence of typical wedge-shaped cystatin features, including the 'QxVxG' and 'Px' motif, as well as the conserved N-terminal glycine residue. Phylogenomic analyses utilizing different cystatin counterparts affirmed the close evolutionary relationship of Bjcystatin-1 and type-1 cystatin homologue. Bjcystatin-1 was predominantly expressed in the gills and hind-gut in a tissue-specific pattern, and its expression was remarkably up-regulated in response to challenge with bacteria or their signature molecules LPS and LTA, suggesting the involvement in immune response. Additionally, the recombinant Bjcystatin-1 (rBjcystatin-1) protein showed significant inhibitory activity towards papain and binding ability to LPS and LTA, indicating its hypothesized role as a pattern recognition receptor in immune response. Subcellular localization results also showed that Bjcystatin-1 was located in the cytoplasm and nucleus, and its overexpression could attenuate the activation of LPS-induced nuclear transcription factors NF-κB. Taken together, our study suggests that amphioxus Bjcystatin-1 acts as a dual role in protease inhibitor and an immunocompetent factor, providing new insights into the immune defense effect of type-1 cystatin in amphioxus.

scRNA-seq analysis of cells comprising the amphioxus notochord.

Cephalochordates occupy a key phylogenetic position for deciphering the origin and evolution of chordates, since they diverged earlier than urochordates and vertebrates. The notochord is the most prominent feature of chordates. The amphioxus notochord features coin-shaped cells bearing myofibrils. Notochord-derived hedgehog signaling contributes to patterning of the dorsal nerve cord, as in vertebrates. However, properties of constituent notochord cells remain unknown at the single-cell level. We examined these properties using Iso-seq analysis, single-cell RNA-seq analysis, and in situ hybridization (ISH). Gene expression profiles broadly categorize notochordal cells into myofibrillar cells and non-myofibrillar cells. Myofibrillar cells occupy most of the central portion of the notochord, and some cells extend the notochordal horn to both sides of the ventral nerve cord. Some notochord myofibrillar genes are not expressed in myotomes, suggesting an occurrence of myofibrillar genes that are preferentially expressed in notochord. On the other hand, non-myofibrillar cells contain dorsal, lateral, and ventral Müller cells, and all three express both hedgehog and Brachyury. This was confirmed by ISH, although expression of hedgehog in ventral Müller cells was minimal. In addition, dorsal Müller cells express neural transmission-related genes, suggesting an interaction with nerve cord. Lateral Müller cells express hedgehog and other signaling-related genes, suggesting an interaction with myotomes positioned lateral to the notochord. Ventral Müller cells also expressed genes for FGF- and EGF-related signaling, which may be associated with development of endoderm, ventral to the notochord. Lateral Müller cells were intermediate between dorsal/ventral Müller cells. Since vertebrate notochord contributes to patterning and differentiation of ectoderm (nerve cord), mesoderm (somite), and endoderm, this investigation provides evidence that an ancestral or original form of vertebrate notochord is present in extant cephalochordates.

Evolution of the gene regulatory network of body axis by enhancer hijacking in amphioxus.

A central goal of evolutionary developmental biology is to decipher the evolutionary pattern of gene regulatory networks (GRNs) that control embryonic development, and the mechanism underlying GRNs evolution. The Nodal signaling that governs the body axes of deuterostomes exhibits a conserved GRN orchestrated principally by Nodal, Gdf1/3, and Lefty. Here we show that this GRN has been rewired in cephalochordate amphioxus. We found that while the amphioxus Gdf1/3 ortholog exhibited nearly no embryonic expression, its duplicate Gdf1/3-like, linked to Lefty, was zygotically expressed in a similar pattern as Lefty. Consistent with this, while Gdf1/3-like mutants showed defects in axial development, Gdf1/3 mutants did not. Further transgenic analyses showed that the intergenic region between Gdf1/3-like and Lefty could drive reporter gene expression as that of the two genes. These results indicated that Gdf1/3-like has taken over the axial development role of Gdf1/3 in amphioxus, possibly through hijacking Lefty enhancers. We finally demonstrated that, to compensate for the loss of maternal Gdf1/3 expression, Nodal has become an indispensable maternal factor in amphioxus and its maternal mutants caused axial defects as Gdf1/3-like mutants. We therefore demonstrated a case that the evolution of GRNs could be triggered by enhancer hijacking events. This pivotal event has allowed the emergence of a new GRN in extant amphioxus, presumably through a stepwise process. In addition, the co-expression of Gdf1/3-like and Lefty achieved by a shared regulatory region may have provided robustness during body axis formation, which provides a selection-based hypothesis for the phenomena called developmental system drift.

Henri de Lacaze-Duthiers and the ascidian hypothesis.

In 1830, Cuvier and Geoffroy Saint-Hilaire confronted each other in a famous debate on the unity of the animal kingdom, which permeated the zoology of the 19th century. From that time, a growing number of naturalists attempted to understand the large-scale relationships among animals. And among all the questions, that of the origin of vertebrates was one of the most controversial. Analytical methods based on comparative anatomy, embryology and paleontology were developed to identify convincing homologies that would reveal a logical sequence of events for the evolution of an invertebrate into the first vertebrate. Within this context, several theories have clashed on the question of the identity of the ancestor of vertebrates. Among the proposals, a group of rather discrete organisms, the ascidians, played a central role. Because he had discovered an ascidian with a particularly atypical larval development, the Molgula, Henri de Lacaze-Duthiers, a rigorous and meticulous naturalist, became involved in the ascidian hypothesis. While the visionary mind of Lacaze-Duthiers led him to establish a particularly innovative methodology and the first marine biology station in Europe, at Roscoff, the tailless tadpole of the Molgula prevented him from recognizing the ancestor of vertebrates. This old 19th century story echoes the ever-present questions driving the field of Eco-Evo-Devo.

Amphioxus as a model to study the evolution of development in chordates.

Cephalochordates and tunicates represent the only two groups of invertebrate chordates, and extant cephalochordates - commonly known as amphioxus or lancelets - are considered the best proxy for the chordate ancestor, from which they split around 520 million years ago. Amphioxus has been an important organism in the fields of zoology and embryology since the 18th century, and the morphological and genomic simplicity of cephalochordates (compared to vertebrates) makes amphioxus an attractive model for studying chordate biology at the cellular and molecular levels. Here we describe the life cycle of amphioxus, and discuss the natural histories and habitats of the different species of amphioxus. We also describe their use as laboratory animal models, and discuss the techniques that have been developed to study different aspects of amphioxus.

Genome-wide transcriptomics and microRNAomics analyses uncover multi-faceted mechanisms to cope with copper stress in ancient macrobenthos amphioxus.

The mechanisms underlying the toxicity of environmental stress are unclear for marine macrobenthos. Copper/Cu has posed the most serious threats to amphioxus, an ancient and model benthic cephalochordate. Herein, a dynamic change in the physiological parameters (GR, SOD, ATP, and MDA) was detected with ROS accumulation in Branchiostoma belcheri exposed to 0.3 mg·L-1 Cu. Transcriptomes and microRNAomes of B. belcheri were generated to investigate the molecular mechanisms by which this amphioxus copes with Cu exposure. Time-specific genes identified at different time points after exposure were involved in the stimulus and immune response, detoxification and ionic homeostasis, aging and the nervous system, sequentially, with prolongation of exposure time, forming a dynamic process of molecular response to Cu stress. In total, 57 differentially expressed miRNAs were identified under Cu stress. Transcriptomics-miRNAomics analyses indicate that these miRNAs targeted genes associated with many key biological processes such as xenobiotics degradation, oxidative stress, and energy metabolism. The constructed miRNA-mRNA-pathway network uncovered a broad post-transcriptional regulatory mechanism in B. belcheri to cope with Cu stress. Overall, this integrated analyses show that enhanced defense response, accelerated ROS elimination, and repressed ATP production constitute a comprehensive strategy to cope with Cu toxicity in the ancient macrobenthos.

Genomic Analysis of Amphioxus Reveals a Wide Range of Fragments Homologous to Viral Sequences.

Amphioxus species are considered living fossils and are important in the evolutionary study of chordates and vertebrates. To explore viral homologous sequences, a high-quality annotated genome of the Beihai amphioxus (Branchiostoma belcheri beihai) was examined using virus sequence queries. In this study, 347 homologous fragments (HFs) of viruses were identified in the genome of B. belcheri beihai, of which most were observed on 21 genome assembly scaffolds. HFs were preferentially located within protein-coding genes, particularly in their CDS regions and promoters. A range of amphioxus genes with a high frequency of HFs is proposed, including histone-related genes that are homologous to the Histone or Histone H2B domains of viruses. Together, this comprehensive analysis of viral HFs provides insights into the neglected role of viral integration in the evolution of amphioxus.

Identification of a novel amphioxus leucine-rich repeat receptor involved in phagocytosis reveals a role for Slit2-N-type LRR in bacterial elimination.

The basal chordate amphioxus is a model for tracing the origin and evolution of vertebrate immunity. To explore the evolution of immunoreceptor signaling pathways, we searched the associated receptors of the amphioxus Branchiostoma belcheri (Bb) homolog of immunoreceptor signaling adaptor protein Grb2. Mass-spectrum analysis of BbGrb2 immunoprecipitates from B. belcheri intestine lysates revealed a folate receptor (FR) domain- and leucine-rich repeat (LRR)-containing protein (FrLRR). Sequence and structural analysis showed that FrLRR is a membrane protein with a predicted curved solenoid structure. The N-terminal Fr domain contains very few folate-binding sites; the following LRR region is a Slit2-type LRR, and a GPI-anchored site was predicted at the C-terminus. RT-PCR analysis showed FrLRR is a transcription-mediated fusion gene of BbFR-like and BbSlit2-N-like genes. Genomic DNA structure analysis implied the B. belcheri FrLRR gene locus and the corresponding locus in Branchiostoma floridae might be generated by exon shuffling of a Slit2-N-like gene into an FR gene. RT-qPCR, immunostaining, and immunoblot results showed that FrLRR was primarily distributed in B. belcheri intestinal tissue. We further demonstrated that FrLRR localized to the cell membrane and lysosomes. Functionally, FrLRR mediated and promoted bacteria-binding and phagocytosis, and FrLRR antibody blocking or Grb2 knockdown inhibited FrLRR-mediated phagocytosis. Interestingly, we found that human Slit2-N (hSlit2-N) also mediated direct bacteria-binding and phagocytosis which was inhibited by Slit2-N antibody blocking or Grb2 knockdown. Together, these results indicate FrLRR and hSlit2-N may function as phagocytotic-receptors to promote phagocytosis through Grb2, implying the Slit2-N-type-LRR-containing proteins play a role in bacterial binding and elimination.

SUMO and PIAS repress NF-κB activation in a basal chordate.

Small ubiquitin-like modifier (SUMO) regulates various biological processes, including the MyD88/TICAMs-IRAKs-TRAF6-NF-κB pathway, one of the core immune pathways. However, its functions are inconsistent between invertebrates and vertebrates and have rarely been investigated in lower chordates, including amphioxus and fishes. Here, we investigated the SUMOylation gene system in the amphioxus, a living basal chordate. We found that amphioxus has a SUMOylation system that has a complete set of genes and preserves several ancestral traits. We proceeded to study their molecular functions using the mammal cell lines. Both amphioxus SUMO1 and SUMO2 were shown to be able to attach to NF-κB Rel and to inhibit NF-κB activation by 50-75% in a dose-dependent fashion. The inhibition by SUMO2 could be further enhanced by the addition of the SUMO E2 ligase UBC9. In comparison, while human SUMO2 inhibited RelA, human SUMO1 slightly activated RelA. We also showed that, similar to human PIAS1-4, amphioxus PIAS could serve as a SUMO E3 ligase and promote its self-SUMOylation. This suggests that amphioxus PIAS is functionally compatible in human cells. Moreover, we showed that amphioxus PIAS is not only able to inhibit NF-κB activation induced by MyD88, TICAM-like, TRAF6 and IRAK4 but also able to suppress NF-κB Rel completely in the presence of SUMO1/2 in a dose-insensitive manner. This suggests that PIAS could effectively block Rel by promoting Rel SUMOylation. In comparison, in humans, only PIAS3, but not PIAS1/2/4, has been reported to promote NF-κB SUMOylation. Taken together, the findings from amphioxus, together with those from mammals and other species, not only offer insights into the functional volatility of the animal SUMO system, but also shed light on its evolutionary transitions from amphioxus to fish, and ultimately to humans.

Three amphioxus reference genomes reveal gene and chromosome evolution of chordates.

The slow-evolving invertebrate amphioxus has an irreplaceable role in advancing our understanding of the vertebrate origin and innovations. Here we resolve the nearly complete chromosomal genomes of three amphioxus species, one of which best recapitulates the 17 chordate ancestor linkage groups. We reconstruct the fusions, retention, or rearrangements between descendants of whole-genome duplications, which gave rise to the extant microchromosomes likely existed in the vertebrate ancestor. Similar to vertebrates, the amphioxus genome gradually establishes its three-dimensional chromatin architecture at the onset of zygotic activation and forms two topologically associated domains at the Hox gene cluster. We find that all three amphioxus species have ZW sex chromosomes with little sequence differentiation, and their putative sex-determining regions are nonhomologous to each other. Our results illuminate the unappreciated interspecific diversity and developmental dynamics of amphioxus genomes and provide high-quality references for understanding the mechanisms of chordate functional genome evolution.

Parallel evolution of amphioxus and vertebrate small-scale gene duplications.

BACKGROUND: Amphioxus are non-vertebrate chordates characterized by a slow morphological and molecular evolution. They share the basic chordate body-plan and genome organization with vertebrates but lack their 2R whole-genome duplications and their developmental complexity. For these reasons, amphioxus are frequently used as an outgroup to study vertebrate genome evolution and Evo-Devo. Aside from whole-genome duplications, genes continuously duplicate on a smaller scale. Small-scale duplicated genes can be found in both amphioxus and vertebrate genomes, while only the vertebrate genomes have duplicated genes product of their 2R whole-genome duplications. Here, we explore the history of small-scale gene duplications in the amphioxus lineage and compare it to small- and large-scale gene duplication history in vertebrates. RESULTS: We present a study of the European amphioxus (Branchiostoma lanceolatum) gene duplications thanks to a new, high-quality genome reference. We find that, despite its overall slow molecular evolution, the amphioxus lineage has had a history of small-scale duplications similar to the one observed in vertebrates. We find parallel gene duplication profiles between amphioxus and vertebrates and conserved functional constraints in gene duplication. Moreover, amphioxus gene duplicates show levels of expression and patterns of functional specialization similar to the ones observed in vertebrate duplicated genes. We also find strong conservation of gene synteny between two distant amphioxus species, B. lanceolatum and B. floridae, with two major chromosomal rearrangements. CONCLUSIONS: In contrast to their slower molecular and morphological evolution, amphioxus' small-scale gene duplication history resembles that of the vertebrate lineage both in quantitative and in functional terms.

Identification of Nodal-dependent enhancer of amphioxus Chordin sufficient to drive gene expression into the chordate dorsal organizer.

The core molecular mechanisms of dorsal organizer formation during gastrulation are highly conserved within the chordate lineage. One of the key characteristics is that Nodal signaling is required for the organizer-specific gene expression. This feature appears to be ancestral, as evidenced by the presence in the most basally divergent chordate amphioxus. To provide a better understanding of the evolution of organizer-specific gene regulation in chordates, we analyzed the cis-regulatory sequence of amphioxus Chordin in the context of the vertebrate embryo. First, we generated stable zebrafish transgenic lines, and by using light-sheet fluorescent microscopy, characterized in detail the expression pattern of GFP driven by the cis-regulatory sequences of amphioxus Chordin. Next, we performed a 5'deletion analysis and identified an enhancer sufficient to drive the expression of the reporter gene into a chordate dorsal organizer. Finally, we found that the identified enhancer element strongly depends on Nodal signaling, which is consistent with the well-established role of this pathway in the regulation of the expression of dorsal organizer-specific genes across chordates. The enhancer identified in our study may represent a suitable simple system to study the interplay of the evolutionarily conserved regulatory mechanisms operating during early chordate development.

A novel nuclear receptor subfamily enlightens the origin of heterodimerization.

BACKGROUND: Nuclear receptors are transcription factors of central importance in human biology and associated diseases. Much of the knowledge related to their major functions, such as ligand and DNA binding or dimerization, derives from functional studies undertaken in classical model animals. It has become evident, however, that a deeper understanding of these molecular functions requires uncovering how these characteristics originated and diversified during evolution, by looking at more species. In particular, the comprehension of how dimerization evolved from ancestral homodimers to a more sophisticated state of heterodimers has been missing, due to a too narrow phylogenetic sampling. Here, we experimentally and phylogenetically define the evolutionary trajectory of nuclear receptor dimerization by analyzing a novel NR7 subgroup, present in various metazoan groups, including cnidarians, annelids, mollusks, sea urchins, and amphioxus, but lost in vertebrates, arthropods, and nematodes. RESULTS: We focused on NR7 of the cephalochordate amphioxus B. lanceolatum. We present a complementary set of functional, structural, and evolutionary analyses that establish that NR7 lies at a pivotal point in the evolutionary trajectory from homodimerizing to heterodimerizing nuclear receptors. The crystal structure of the NR7 ligand-binding domain suggests that the isolated domain is not capable of dimerizing with the ubiquitous dimerization partner RXR. In contrast, the full-length NR7 dimerizes with RXR in a DNA-dependent manner and acts as a constitutively active receptor. The phylogenetic and sequence analyses position NR7 at a pivotal point, just between the basal class I nuclear receptors that form monomers or homodimers on DNA and the derived class II nuclear receptors that exhibit the classical DNA-independent RXR heterodimers. CONCLUSIONS: Our data suggest that NR7 represents the "missing link" in the transition between class I and class II nuclear receptors and that the DNA independency of heterodimer formation is a feature that was acquired during evolution. Our studies define a novel paradigm of nuclear receptor dimerization that evolved from DNA-dependent to DNA-independent requirements. This new concept emphasizes the importance of DNA in the dimerization of nuclear receptors, such as the glucocorticoid receptor and other members of this pharmacologically important oxosteroid receptor subfamily. Our studies further underline the importance of studying emerging model organisms for supporting cutting-edge research.

Localization of trypsin, chymotrypsin and elastase in the digestive tract of amphioxus Branchiostoma japonicum with implications to the origin of vertebrate pancreas.

The organ pancreas is characteristic of all vertebrates, but its evolutionary origin remains largely unknown. The serine proteinases trypsin, chymotrypsin and elastase are produced primarily in the pancreas in vertebrates. Here we clearly demonstrate by enzymatic activity assay, histochemical staining and in situ hybridization that trypsin, chymotrypsin and elastase are widely distributed in the digestive tract of the amphioxus Branchiostoma japonicum, especially in the hepatic caecum and mid-gut. This suggests that an extensive region of the amphioxus digestive tract including the hepatic caecum is homologous to the vertebrate exocrine pancreatic cells, providing a new angle for the study of the origin and evolution of vertebrate pancreas.

Three-dimensional fine structure of fibroblasts and other mesodermally derived tissues in the dermis of adults of the Bahamas lancelet (Chordata, Cephalohordata), as seen by serial block-face scanning electron microscopy.

Tissues of adult cephalochordates include sparsely distributed fibroblasts. Previous work on these cells has left unsettled such questions as their developmental origin, range of functions, and even their overall shape. Here, we describe fibroblasts of a cephalochordate, the Bahamas lancelet, Asymmetron lucayanum, by serial block-face scanning electron microscopy to demonstrate their three-dimensional (3D) distribution and fine structure in a 0.56-mm length of the tail. The technique reveals in detail their position, abundance, and morphology. In the region studied, we found only 20 fibroblasts, well separated from one another. Each was strikingly stellate with long cytoplasmic processes rather similar to those of a vertebrate telocyte, a possibly fortuitous resemblance that is considered in the discussion section. In the cephalochordate dermis, the fibroblasts were never linked with one another, although they occasionally formed close associations of unknown significance with other cell types. The fibroblasts, in spite of their name, showed no signs of directly synthesizing fibrillar collagen. Instead, they appeared to be involved in the production of nonfibrous components of the extracellular matrix-both by the release of coarsely granular dense material and by secretion of more finely granular material by the local breakdown of their cytoplasmic processes. For context, the 3D structures of two other mesoderm-derived tissues (the midline mesoderm and the posteriormost somite) are also described for the region studied.

Highly distinct genetic programs for peripheral nervous system formation in chordates.

BACKGROUND: Vertebrates develop their peripheral nervous system (PNS) from transient unique embryonic structures, the neural crest, and the ectodermal placodes that are located at the border of the forming central nervous system. By contrast, in the invertebrate chordates, amphioxus and ascidians, a large part of the PNS originates at the opposite of the embryo, in the ventral ectoderm. In both groups, a biphasic mechanism regulates ventral PNS formation: high BMP levels specify a neurogenic territory within which glutamatergic epidermal sensory neuron formation is controlled by the Notch pathway. Given these similarities and the phylogenetic relationships within chordates, it is likely that ventral PNS is an ancestral feature in chordates and that it has been lost in vertebrates. RESULTS: In order to get insights into the molecular control of ventral PNS formation and to test the hypothesis of their homology and potential contribution to the emergence of vertebrate PNS, we undertook a close comparison of ventral PNS formation in the ascidian Phallusia mammillata and the amphioxus Branchiostoma lanceolatum. Using timed RNA-seq series, we identified novel markers of the ventral PNS during different phases of its development in both species. By extensively determining the expression of paralogous and orthologous genes, we observed that only a minority of genes have a shared expression in the ventral PNS. However, a large fraction of ventral PNS orthologous genes are expressed in the dorsally forming PNS of vertebrates. CONCLUSIONS: Our work has significantly increased the molecular characterization of ventral PNS formation in invertebrate chordates. The low observed conservation of gene expression in the ventral PNS suggests that the amphioxus and ascidian ventral PNS are either not homologous, or alternatively extensive drift has occurred in their regulatory mechanisms following a long period (600 My) of separate evolution and accelerated evolution in the ascidian lineage. The homology to genes expressed in the dorsally forming PNS of vertebrates suggests that ancestral sensory neurons gene networks have been redeployed in vertebrates.

Joint profiling of gene expression and chromatin accessibility during amphioxus development at single-cell resolution.

Vertebrate evolution was accompanied by two rounds of whole-genome duplication followed by functional divergence in terms of regulatory circuits and gene expression patterns. As a basal and slow-evolving chordate species, amphioxus is an ideal paradigm for exploring the origin and evolution of vertebrates. Single-cell sequencing has been widely used to construct the developmental cell atlas of several representative species of vertebrates (human, mouse, zebrafish, and frog) and tunicates (sea squirts). Here, we perform single-nucleus RNA sequencing (snRNA-seq) and single-cell assay for transposase accessible chromatin sequencing (scATAC-seq) for different stages of amphioxus (covering embryogenesis and adult tissues). With the datasets generated, we constructed a developmental tree for amphioxus cell fate commitment and lineage specification and characterize the underlying key regulators and genetic regulatory networks. The data are publicly available on the online platform AmphioxusAtlas.

PACAP/GCGa Is an Important Modulator of the Amphioxus CNS-Hatschek's Pit Axis, the Homolog of the Vertebrate Hypothalamic-Pituitary Axis in the Basal Chordates.

The Hatschek's pit in the cephalochordate amphioxus, an invertebrate deuterostome basal to chordates is suggested to be the functional homolog structure of the vertebrate adenohypophysis based on anatomy and expression of homologous neuroendocrine genes. However, the endocrine potential of the cephalochordate Hatschek's pit remains to be demonstrated as well as the physiological actions of the secreted neuropeptides. In this study, we have explored the distribution and characterize the potential function of the amphioxus PACAP/GCG precursor, which is the ortholog of the hypothalamic PACAP neuropeptide in vertebrates. In amphioxi, two PACAP/GCG transcripts PACAP/GCGa and PACAP/GCGbc that are alternative isoforms of a single gene with different peptide coding potentials were isolated. Immunofluorescence staining detected their expression around the nucleus of Rohde, supporting that this structure may be homologous of the neurosecretory cells of the vertebrate hypothalamus where abundant PACAP is found. PACAP/GCGa was also detected in the infundibulum-like downgrowth approaching the Hatschek's pit, indicating diffusion of PACAP/GCGa from the CNS to the pit via the infundibulum-like downgrowth. Under a high salinity challenge, PACAP/GCGa was upregulated in amphioxi head and PACAP/GCGa treatment increased expression of GHl in Hatschek's pit in a dose-dependent manner, suggesting that PACAP/GCGa may be involved in the regulation of GHl via hypothalamic-pituitary (HP)-like axis similar as in the vertebrates. Our results support that the amphioxus Hatschek's pit is likely to be the functional homolog of pituitary gland in vertebrates.